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Abstract 

To  follow  the  CCTF  2006  recommendation  of  using  high-precision  GPS  carrier-phase 
data  (CP)  for  TAI  generation,  BIPM  computes  monthly  the  GPS  PPP  (Precise  Point 
Positioning)  time  links  and  publishes  the  result  on  the  TAI  ftp  site.  The  GPS  PPP  is  a 
combined  solution  of  the  codes  and  CP  information.  The  latter  assigns  the  high  precision  to 
the  PPP  result.  It  is  shown  that  the  PPP  links  have  a  better  time  stability  than  the  TW  ones 
in  short  and  middle  terms  (up  to  a  week).  On  the  other  hand,  TW  could  be  more  stable  on 
the  long  term,  due  to  its  symmetric  measurement  procedure  and  its  more  accurate  absolute 
calibration.  The  combination  of  PPP  and  TW  may  benefit  the  advantages  of  both  links. 
Because  the  GPS  and  TW  are  completely  independent  techniques,  combination  of  them 
makes  the  results  less  noisy,  less  biased,  and  more  robust.  The  combination  gives  also  a 
possibility  to  analyze  the  errors  presented  in  each  technique.  One  of  them  is  the  diurnal 
deviation  in  some  TW  KU  band  links.  In  this  paper,  we  first  briefly  review  the  recent  study 
on  the  uncertainties  of  PPP  and  TW,  then  quickly  recall  the  model  used  for  the  combination, 
and  finally  evaluate  of  the  method  proposed  using  measuring  data.  The  TAI  time  network  is 
highly  redundant.  This  approach  gives  a  new  idea  about  fully  using  the  potential  of  the  two 
fundamental  time  transfer  techniques:  GPS  and  TW. 


1.  INTRODUCTION 

GPS  PPP  stands  for  time  transfer  using  the  Precise  Point  Positioning  technique.  Many  authors  studied 
this  technique  and  concerned  problems,  such  as  Kouba,  Heroux,  Ray,  Cerretto,  Orgiazzi,  Guyennon, 
Dach,  Deffaigne,  Petit  and  Jiang,  etc.  [1-5].  PPP  time  transfer  is  now  becoming  operational.  BIPM 
gathered  the  study  results  and  reported  to  the  CCTF  2006.  It  is  recognized  that  the  GPS  carrier-phase 
information  is  a  powerful  tool  to  improve  the  quality  of  time  links  and  to  answer  the  new  challenge  of 
comparing  new  generation  clocks.  Following  accordance  with  the  CCTF  2006  Recommendation  4, 
BIPM  will  routinely  compute  the  PPP  time  transfer  and  put  available  the  results  on  its  ftp  site. 

We  know  that  the  stability  in  the  PPP  solution  is  dominated  by  the  accuracy  and  high  rate  of  the  carrier- 
phase  information.  Earlier  studies  have  proven  that  the  combination  of  TW  and  the  scale-less  GPS  carrier 
phase  can  considerably  improve  the  TW  time  transfer  stability  without  changing  its  calibration  [4,5].  In 
this  paper,  we  will  show  that  PPP  may  be  used  in  the  same  way  with  similar  result  for  time  and  frequency 
transfer. 

In  the  TAI  computation,  about  20%  of  TAI  laboratories  which  concentrate  nearly  70%  of  the  total  weight 
of  the  TAI  clocks  and  host  all  1 1  PFS  (Primary  Frequency  Standards),  have  both  TW  and  PPP  data.  In 
fact,  TAI  time  transfer  network  is  highly  redundant.  Almost  all  the  TW  laboratories  are  backed  up  by  the 
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Public  reporting  burden  for  the  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources,  gathering  and 
maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of  information, 
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geodetic  GPS  receivers.  However,  the  current  TAI  time  link  strategy  is  that  we  use  only  the  primary 
technique,  either  TW  or  GPS,  but  never  both. 

In  view  of  time  stability,  TW  and  PPP  show  their  advantages  in  different  averaging  times.  We  show  this 
by  comparing  the  uncertainty  types  A  and  B  (uA,  uB)  in  TAI  time  transfer.  First  uA  represents  the 
statistical  uncertainty  and  is  dominated  by  the  measurement  error.  For  short  terms  up  to  several-day 
averaging,  the  stability  of  PPP  is  obviously  better  than  that  of  TW  due  to,  as  discussed  below,  the  rather 
bigger  noises  and  diumals  in  TW  measurements.  For  middle  terms  up  to  month(s),  it  is  difficult  to 
discriminate  both  techniques,  but  TW  might  be  more  advantageous  than  any  GPS  techniques.  Secondly, 
uB  represents  the  uncertainty  of  the  absolute  calibration  of  the  equipment.  The  uB  of  PPP  is  that  of  the 
P3  code  and  is  estimated  to  be  of  order  5  ns  for  a  link.  The  uB  of  TW  may  be  as  small  as  1  ns.  However, 
the  TW  calibration  also  depends  on  the  communication  satellites  involved  and  any  change  in  the  ground 
equipment  or  in  the  satellite  would  lose  the  TW  calibration.  Obviously,  each  technique  has  its  advantages 
and  disadvantages.  Combining  the  two  techniques  makes  sense:  a)  They  involve  completely  independent 
hardware  and  use  different  satellite  systems;  b)  The  combination  allows  one  to  keep  the  advantages,  to 
reduce  the  disadvantages,  and  to  make  the  result  less  noisy,  less  biased,  and  more  robust;  c)  The 
combination  gives  a  possibility  of  recovering  the  TW  calibration  at  the  PPP  resolution  of  about  100  ps; 
that  is,  we  can  first  calibrate  the  GPS  receiver  with  respect  to  the  TW  and  then  when  the  TW  calibration  is 
lost,  we  can  restore  it  with  the  GPS  (this  is  the  so-called  bridge  calibration),  and  d)  It  is  a  solution  for  the 
high  redundancy  in  the  TAI  network  so  as  to  fully  use  the  total  potential  of  different  time  link  data. 
Furthermore,  we  cannot  evade  this  problem  as  more  and  more  new  techniques  are  coming,  e.g.,  Glonass, 
Galileo,  Compass,  etc.;  the  question  of  how  to  combine  different  techniques  is  inevitable. 

The  mathematical  method  used  is  the  Vondrak-Cepek  combined  smoothing  [7],  which  has  been  discussed 
in  the  earlier  publication  [4,5].  We  review  here  only  the  principle  of  the  method.  The  input  required  is  of 
two  series:  SI  is  the  TW  measurements,  which  are  the  clock  differences,  and  S2  is  the  derivatives  of  the 
clock  differences  calculated  with  the  rate  of  PPP.  By  this  configuration  of  the  data  sets,  we  know,  onto 
the  output,  the  TW  will  contribute  to  the  calibration  with  its  uncertainty  uB  and  the  GPS  PPP  will 
contribute  mainly  the  improvement  in  the  time  stability  with  its  uncertainty  uA.  The  combined 
smoothing  is  performed  under  the  three  conditions:  F:  close  to  SI;  F°:  close  to  S2,  and  S:  smoothness  of 
the  combined  series  S3,  written  as: 

Q  =  S  +  eF  +  s°F°  =  min 

Here  s  and  s°  are  the  smoothing  coefficients.  By  adjusting  them,  we  obtain  the  optimal  combination 
result.  It  is  an  iteration  procedure  based  on  the  spectral  analysis  of  the  TW  and  PPP  data  sets. 


2.  NUMERICAL  TESTS 

We  use  two  tools  to  examine  the  gains  obtained  by  the  TW-PPP  combination.  The  first  is  the  well  known 
Time  deviation  Tdev,  which  measures  the  stability,  and  the  second  is  the  closure  analysis,  which  also 
addresses  (to  some  level)  the  accuracy.  We  first  discuss  the  closure  analysis  as  an  index  of  the  time  link 
quality  and  then  discuss  in  detail  how  much  the  gains  will  be.  Because  the  TW  calibration  is  not  a  point 
in  the  following  tests,  we  use  all  the  available  TW  data,  whether  they  are  calibrated  or  not.  Three 
triangles,  all  with  H-maser  links,  are  fully  analyzed.  They  are  the  short  and  long  baselines  between  PTB 
in  Germany,  IT  in  Italy,  NIST  and  USNO  in  USA,  and  the  triangle  between  KRIS  in  South  Korea,  TL  in 
Chinese  Taipei,  and  NICT  in  Japan.  TW  and  GPS  PPP  data  were  gathered  as  a  typical  TAI  month:  TAI 
0709  (Sept.  2007,  MJD  54344-54369,  totally  25  days). 
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Table  1.  Statistics  of  the  TW0709  closures. 
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2.1  The  TW  Link  Closure 

The  definition  of  the  TW  link  closure  is  simple:  a  triangle  closure  equals  the  sum  of  the  three  time  link 
vectors  between  the  three  time  laboratories  (Fig.  1): 


Fig.  1.  Definition  of  a  triangle  closure  of  TW  time  links: 

Closure  =  TW1+TW2+TW3 

=  (Clock2-Clock  1  )+(Clock3-Clock2  )+(Clock  1  -Clock3 ) 
=  V1+V2+V3 


VI,  V2,  and  V3  are  the  measurement  errors  of  the  time  links  TW1,  TW2,  and  TW3  measured 
independently.  By  definition,  all  the  clocks  are  cancelled  in  the  closure,  if  all  measurements  are 
simultaneous.  Because  the  clocks  are  H-masers,  any  effect  of  nonsimultaneity  is  negligible  and  the 
closure  should  equal  to  zero  if  there  were  not  measurement  errors.  A  nonzero  closure  is  a  true  error.  The 
closures  statistic  gives  us  a  true  error  index  to  judge  the  quality  of  the  TW  links.  On  average,  there  are 
10-12  points  per  day  and  such  about  260  triangle  closures  for  the  TAI0709.  In  the  above  equation,  the 
terms  VI,  V2,  and  V3  can  be  rigorously  and  uniquely  determined.  But  this  is  not  the  point  of  this  paper. 
Interested  people  can  refer  to  [8]. 

Table  1  displays  the  statistics  of  the  55  TW  closures  for  TAI  0709.  The  mean  values  are  not  zero,  but 
vary  between  ±  2  ns.  The  STD  (standard  deviation  around  the  mean  value)  of  the  closures  varies  between 
±  1  ns,  except  for  a  few  cases. 
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Mean=  -909  ps 
RMS=  974  ps 
STD=  348  ps 


+  +  +  + 
+  +  +  + 
+  +  +  + 


+  +  +  + 


+  +  +  + 


+  +  +  + 
+  +  +  + 
+  +  +  + 


OOOOOOOOOO 

fflOOOOOOOOO 

Wcor-'X'Ln^rcocMHO 


DOOOOOOOOO 
300000000 
or-OLn'TrocsjH 


Fig.  2.  Histogram  of  the  TW  closures  of 
the  TW  link  triangle  NIST-IT-PTB.  It  is 
of  the  Europe-America  long-distance 
baseline.  There  are  10-12  points  per  day 
and  totally  236  points  over  the  25  days 
for  TAI  0709,  2007  Sept. 


Fig.  2  is  the  histogram  of  the  TW  closures  for  the  baseline  triangle  NIST-IT-PTB.  The  mean  value  is 
-909  ps  with  the  RMS  947  ps.  The  STD  is  348  ps  due  to  mainly  the  measurement  noises  and  diurnal 
errors. 

2.2  Measurement  Noises  and  Diurnals 


< -  25  days  - > 


Fig.  3.  Comparison  of  the  meas¬ 
urement  noises  between  the  TW 
and  PPP  for  the  TAI  0709  links  IT- 
NIST  and  IT-USNO,  all  with  HM 
clocks.  The  peak-to-peak  scatter  is 
2-3  ns  for  TW  and  0.3  ns  for  PPP. 
It  might  not  be  a  normal  case  for 
the  TW  links,  cf.  Fig.  4. 


Fig.  3  is  the  comparison  of  the  noise  level  between  TW  and  PPP  time  links  IT-NIST  and  IT-USNO.  They 
are  both  the  long-distance  and  Fl-maser  links.  The  black  points  are  of  TW  with  a  peak-to-peak  scatters  of 
2-3  ns,  which  are  obviously  due  to  noise  in  this  link.  Tdev  analysis  (Fig.  4)  shows  the  level  of  diurnal 
signals,  which  dominate  up  to  averaging  time  of  about  2  days.  The  blue  points  are  of  the  PPP,  with  a 
scatter  about  several  hundred  ps. 

Several  authors  have  already  discussed  the  diurnals.  At  present,  we  cannot  give  a  reasonable  explanation 
for  the  physical  mechanism  of  the  diurnals.  Therefore,  we  have  not  yet  a  successful  method  to  remove  or 
reduce  the  diurnals  in  the  raw  TW  data.  What  we  can  do  is  to  remove  or  reduce  them  with  the  help  of 
GPS  PPP.  The  curves  on  the  left  in  Fig.  4  show  that  strong  diurnals  are  clearly  visible  throughout  the 
whole  month.  Further  proof  can  be  seen  in  the  Tdev.  In  the  middle  plots,  one  cannot  miss  the  strong 
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diurnal  signals  in  the  Tdev.  On  the  other  hand,  the  plots  on  the  right  are  the  Tdev  of  the  related  PPP  links, 
with  only  very  slight  diumals  that  can  so  far  be  ignored. 


Raw  TW  links  2007 

< -  25  days  - 
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Fig.  4.  The  diumals  in  TW  time  links.  Left:  the  raw  TW  links  where  clear  diumals  are 
visible;  Middle:  the  Tdev  corresponding  to  the  raw  TW  links;  Right:  the  raw  PPP  links, 
where  no  significant  diumals  exist.  Units  are  ns  for  the  links  and  0. 1  ns  for  Tdev. 


2.3  Combination  of  TW  and  PPP 

By  the  above  comparisons  and  in  view  of  the  noise  and  diumals,  the  PPP  is  more  advantageous  than  TW. 
This  is  the  physical  basis  of  the  combination  of  the  two  techniques. 


TW+PPP  links  Combined/ns 
< -  25  days  - > 


TDev  of  the  combined 

dq 


0/ 


0;0  h/2  h  d/8  d/4  d/2  day  ddd  Vk 


h/2  h  d/8  d/4  d/2  djy 


Fig.  5.  Results  of  the  combined  TW+PPP  links.  Units  are  ns  for  the  links  and  0. 1  ns  for  Tdev. 


Fig.  5  is  the  combined  time  link.  Compared  the  curves  in  Fig.  5  with  the  same  TW  links  and  the  related 
Tdev  in  Fig.  4,  both  the  noises  and  the  diumals  are  greatly  reduced. 
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Now  we  look  at  the  combination  using  closures  as  another  tool.  As  mentioned  above,  the  nonzero  closure 
represents  a  true  error.  The  closure  analysis  gives  a  true  index  of  how  much  the  improvement  in  the 
combined  result.  Fig.  6  presents  the  histograms  of  the  closures  of  the  three  triangles.  The  blue 
histograms  are  of  the  raw  TW  closures  and  the  red  ones  are  of  the  closures  of  the  combined  links.  The 
mean  value  and  the  RMS  have  practically  no  change.  This  implies  the  combination  keeps  the  calibration 
defined  by  TW.  The  STD  is  reduced  about  40%  for  NIST-IT-PTB  and  as  big  as  476%  for  USNO-IT- 
PTB.  The  Asia  triangle  NICT-KRIS-TL  is  relatively  smaller.  The  STD  of  the  raw  data  closures  is 
already  very  small;  only  121  ps  that  approximate  the  limit  of  the  TW  stability.  But  even  in  this  case,  the 
STD  is  reduced  34%.  This  further  proves  that  the  time  stability  of  PPP  is  better  than  that  of  the  best  TW 
links.  The  closures,  especially  those  with  large  standard  deviation,  are  greatly  reduced.  Comparing  Fig. 
4,  5,  and  6,  closures  are  reduced  partially  because  the  diumals  become  smaller. 


Closure  A:  NICT-KRIS-TL 


Raw : 

Combined  : 

Mean=  -577  ps 

Mean=  -570  ps 

RMS=  589  ps 

RMS=  577  ps 

STD=  121  ps 

STD=  90  ps 

STD  reduced  34% 


Fig.  6.  Comparison  between  the  raw 
TW  closures  and  the  TW-PPP 
combined  closures.  Units  are  ps.  The 
mean  and  RMS  keep  no  change 
statistically.  This  implies  the  TW 
calibrations  (uB)  have  no  change  after 
the  combination.  The  STD  is  greatly 
reduced.  The  quality  of  the  time  link  is 
considerately  improved  in  time  stability 
(uA). 
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2.4  The  Redundancy  and  the  Combination 

Since  last  Oct  2006,  the  TAI  time  transfer  network  has  become  a  one-pivot  system;  the  pivot  is  PTB. 
Flowever,  there  are  many  TW  laboratories  that  have  no  direct  TW  link  with  PTB.  Their  TW  results 
cannot  be  used  for  TAI.  An  example  is  the  Asian  lab  TL.  To  compute  the  TW  link  TL-PTB,  we  have  to 
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use  the  NICT  as  an  intermediate  lab  to  bridge  the  time  transfer  (Fig.  7).  And  on  the  other  hand,  all  PPP 
links  are  available  and  their  quality  is  independent  of  the  distance.  An  application  of  the  proposed 
method  is  to  combine  two  TW  links  with  a  PPP  link,  which  comprise  a  triangle.  The  comparison  of  the 
Tdev  of  the  TW  and  PPP  in  Fig.  7  shows  that  there  is  no  doubt  this  combination  will  considerably 
improve  the  quality  of  this  Asia-Europe  long-distance  link.  This  is  only  an  example  of  the  redundancy  of 
the  TAI  time  transfer  data.  In  fact,  on  many  baselines,  there  are  300%  redundant  data  in  the  TAI  network. 
Combination  of  TW  and  PPP  is  a  solution  that  fully  uses  the  redundancy. 


NICT 


Fig.  7.  TL-PTB  =  (TL-NICT)  -  (PTB-NICT):  Combining  two  indirect  TW  links  and  a 
direct  PPP  link.  The  Tdev  of  the  PPP  is  well  below  that  of  the  TW.  The  combination 
should  greatly  improve  the  quality  of  the  Asia-Europe  long-baseline  time  transfer  and  is  a 
solution  of  high  redundancy  in  the  TAI  network. 

3.  CONCLUSION 

Combining  two  independent  links  TW  +  PPP  allows:  1)  full  use  of  the  advantages  of  the  calibration  and 
the  long-term  stability  in  TW  and  the  short-term  stability  in  PPP  and  2)  reduction  of  the  short-term  noises 
and  the  diumals  in  TW  as  well  as  the  poor  calibration  and  possible  biases  in  PPP.  The  combined  result  is 
of  higher  short-  and  long-term  stability  and  more  robust. 

20%  of  TAI  laboratories  have  both  TW  and  PPP.  They  transfer  2/3  of  the  TAI  clocks  and  all  the  1 1  PFS. 
Combining  TW+PPP  of  these  labs  will  effectively  improve  TAI. 


The  proposed  method  allows  using  the  high  redundancy  in  the  worldwide  TAI  TW-GPS  time  transfer 
network. 
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